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Abstract 

GUT monopoles captured by the Sun's gravitation are expected to catalyze proton 
decays via the Callan-Rubakov process. In this scenario, protons, which initially 
decay into pions, will ultimately produce fg, i'^ and P^. After undergoing neutrino 
oscillation, all neutrino species appear when they arrive at the Earth, and can be 
detected by a 50,000 metric ton water Cherenkov detector, Super-Kamiokande (SK). 
A search for low energy neutrinos in the electron total energy range from 19 to 55 
MeV was carried out with SK and gives a monopole flux limit of FM((7o/lmb) < 
6.3 X 10-2''(/3m/10~3)2 cm~2 g-^ sr'^ at 90% C.L., where /3m is the monopole 
velocity in units of the speed of light and ctq is the catalysis cross section at /3m = 1- 
The obtained limit is more than eight orders of magnitude more stringent than the 
current best cosmic-ray super massive monopole flux limit, Fm < 1 x 10~^^ cm~^ 
s~^ sr~^ for /3m < 10~^ and also two orders of magnitude lower than the result of 
the Kamiokande experiment, which used a similar detection method. 

Key words: GUT monopole, neutrino, water Cherenkov detector 
PACS: 14.80.Hv, 25.30.Pt, 29.40.Ka 



1 Introduction 



Grand Unified Theories (GUTs) predict superheavy magnetic monopoles (GUT 
monopoles) produced in the very early universe [T||2] . GUT monopoles are 
predicted to have appeared as topological defects at the phase transition of 
vacuum, where the GUT gauge group spontaneously broke to leave the U(l) 
of electromagnetism. If we take as an example the temperature of the phase 
transition at which the monopoles were formed to equal 10^^ GeV and as- 
sume the average production rate of about one monopole per horizon at that 
time [3] , the density of the monopoles today would exceed the critical density 
of the universe by more than 14 orders of magnitude. Even if we circumvent 
this problem by resorting to inflationary universe scenarios [HE], we cannot 
avoid a large uncertainty on the monopole flux in the universe since the flux 
depends on such parameters as monopole mass and the reheating tempera- 
ture. In fact, due to the wide variety of elementary particle models, several 
models are compatible with the level of the Parker bound (~ 10~^^ cm~^ 
g-i g];-i) [5|7|5] . and a flux in that range can be relatively easily detected 
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by underground experiments. Arafune et al. [9] pointed out that copious low 
energy neutrinos might be emitted when monopoles accumulating inside the 
Sun catalyze proton decays, 

p ^ (p°, w, r^, ir+, ■ ■ ■ ) + e+(or /i+) (1) 



along their paths with cross sections typical of strong interactions via the 
Callan-Rubakov process (TOUT]. When decay mesons produced by the above 
process subsequently d ecay into positive pions, ((p°, cj, r^, K^ , ■■■)—!• 7r+), z/g, 



.„ and P„ aie produce<inIby the reactions 



7r+ -^ /i+ + V^ 

/i+ -> e+ + z/e + v^ 



(2) 
(3) 
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Fig. 1. Expected neutrino spectra at the SK detector. Neutrino oscillation effects 
are taken into account, /i, r components are collectively denoted as h since these two 
flavors of neutrinos have identical cross sections for electron scattering. The spikes 
at 29.79 MeV originate from u^ produced by the two-body decay of vr^ in Eq. (2), 
which also oscillates into V(.- 

After undergoing neutrino oscillation, all neutrino species are present when 
they arrive at the Earth (Fig. [T]), and such neutrinos can be detected by a 
water Cherenkov detector. Throughout this paper, the neutrinos and antineu- 
trinos are assumed to oscillate with the same parameters of sin^ 6'i2=0.31, 



^^ Under the assumption that monopoles are accumulated in the center of the Sun, 
pions decaying in flight are negligible due to the high density of hydrogen in the 
Sun. 



sin^ ^13=0.02 and A?7i2i=7.6 x 10 ^ eV^ [12]. The uncertainties of the mixing 
angles are taken into account in the systematic errors. 



2 Detector and Expected signal 



We searched for monopole- induced neutrinos using Super-Kamiokande (SK) [T3] . 
a large water Cherenkov detector located inside Japan's Kamioka mine. SK is 
a high-performance neutrino detector consisting of 50,000 metric tons of pure 
water. Since there are 7-ray backgrounds near the wall, the fiducial volume 
for this search, which amounts to 22,500 metric tons, is defined to be more 
than two meters from the walls of the inner detector (ID). 

Monopole-induced neutrinos include all six types, and so for the monopole- 
induced neutrino search both electron elastic scattering, i'r^{i?x)+e~ — )■ i'x{i^x) + 
e~ , and inverse beta decay, i^e +P — ^ e^ + n, were assumed to contribute. Fig- 
ure [2] shows the expected spectra of elastic scattering recoil electrons and in- 
verse beta decay positrons in SK. The Mikheyev-Smirnov-Wolfenstein (MSW) 
effects [13] on oscillation in both the Sun and the Earth were also taken 
into account. We adopted a simple approximate analytic formula presented 
in Ref. [15] for the calculation of the MSW effects and assumed the adiabatic- 
ity condition is fulfilled, i.e., the density changes so slowly inside the Sun and 
the Earth that the neutrino mass eigenstates propagate independently. The 
density profiles of the Sun and the Earth were taken from Refs. fIB] and [T7] , 
respectively. 



3 Data reduction 



In this analysis, we used 2853 live days of data consisting of SK-1 (Apr. 1996 - 
Jul. 2001 : 1497.4 days), SK-11 (Oct. 2002 - Oct. 2005 : 793.7 days) and SK-111 
(Jul. 2006 - Aug. 2008 : 562.0 days). In SK-11 the energy resolution was worse 
because the number of photomultiplier tubes (PMTs) was about half as many 
as those in the other phases. This difference in detector geometry was taken 
into account in the calculation of the expected spectra. However, as previously 
demonstrated by the SK solar neutrino analysis [18j, the broad energy spectra 
typical of monopole-induced neutrinos are hardly affected by the resolution. 

The background events of the monopole-induced neutrino events are mainly 
caused by the atmospheric neutrinos, the solar neutrinos and muon-induced 
spallation products. To remove the spallation products and the solar neutrino 
events, we set the lower energy thresholds to 19 MeV for SK-1 and 111, and 



20 MeV for SK-II in this analysia ^^ I. The upper energy threshold was set to 
55 MeV, which was determined by the end-point of the recoil electron 
positron) spectrum (Fig. [2]). 
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Fig. 2. Expected energy spectra of elastic scattering recoil electron (left) and inverse 
beta decay positron (right) in SK assuming neutrino flux of 3.0 x 10^ cm~^s~^ 
without energy resolution included. The discontinuity of some spectra at 29.79 MeV 
is caused by the spikes shown in the expected neutrino energy spectra in Fig. [TJ h 
represents n and r components. 

Besides setting the fiducial volume and energy criteria, some additional back- 
ground cuts were applied to the data. After the event selection that removes 
cosmic ray muons and detector noise events, the data sample is subjected to a 
spallation cut. Cosmic-ray muons can spall oxygen and create unstable nuclei 
called spallation products (/i + ^^O — )■ /x + X). This is one of the most abun- 
dant backgrounds in the < 20 MeV region, and as mentioned above, the ability 
to remove the background events mostly determines the lower threshold for 
the monopole-induced neutrino search. The spallation background events are 
reduced by a likelihood method that uses timing, position and photo-electron 
information of the muons preceding the candidate events. 

Then, we applied the Cherenkov opening angle cut to the remaining data sam- 
ple. The opening angle is estimated by histogramming all the angles uniquely 
determined by the reconstructed vertex and the possible combinations of three 
hit PMTs. The histogram is divided into 100 angle bins and the peak is lo- 
cated by finding the successive seven bins with the largest number of entries. 
The angle corresponding to the midpoint of the seven bins is regarded as the 
Cherenkov angle of the event. Most of the remaining visible atmospheric (anti- 
)muon neutrino events are removed by this cut. Electrons with E > 18 MeV 



T^ In the SK solar neutrino analysis, we reconstruct the event energy assuming all 
Cherenkov photons in an event come from a single electron, and we define "energy" 
as the total energy (not the kinetic energy) of the electron. Throughout this paper, 
we use this definition. 



have a Cherenkov angle 9c of about 42°, while visible muons remaining in 
the signal energy range have momenta less than ~ 250 MeV/c, which corre- 
sponds to the maximum 9c value of 36°. Thus, taking into account the finite 
resolution of the Cherenkov angle measurements, events with 9c < 38° were 
removed. Also we removed events with 9c > 50°. Those events do not have 
a clear Cherenkov ring pattern and originate from multiple 7-rays emitted 
by excited nuclei created in neutral current (NC) interactions of atmospheric 
neutrinos. 

Some events originating from the outside of the fiducial volume are recon- 
structed within the fiducial volume of SK. These events are 7-rays from the 
materials of the detector structure and the surrounding rock. To remove such 
events, we cut on the event's distance to the ID wall projected backwards from 
its vertex position along its reconstructed direction. Events with this distance 
less than 300 cm were removed (7-ray cut [18j) 



t3 



We removed charged tt background events based on the sharpness of the 
Cherenkov ring pattern. As a charged vr interacts with or is absorbed by 
a nucleus and is less influenced by multiple Coulomb scattering than elec- 
trons, it generates a sharp Cherenkov ring pattern. Events having two or more 
Cherenkov rings were also removed. 

The remaining backgrounds include low energy muons produced by atmo- 
spheric z/^ charged current (CC) interactions and their decay electrons. Those 
events are removed using time and spatial correlations; any event pairs having 
a time difference and a spatial distance less than 50 fis and 500 cm, respec- 
tively, were removed. In case the candidate event is the second one in a pair, we 
removed event pairs with the first events having total photo-electrons (p.e.) 
more than 1000 p.e. for SK-I and III, and 500 p.e. for SK-II, which corre- 
spond to 90 ~ 100 MeV electron equivalent energy. Although this selection 
(decay electron cut) removes some fraction of the atmospheric neutrino events, 
there still remain decay electrons which are not accompanied by nuclear de- 
excitation 7-rays, and whose parent muons are invisible because their energies 
are below Cherenkov threshold. Those background events are distinguished 
from the signal events using angular correlation with the Sun direction, which 
will be described later. 



4 Analysis and Results 



The numbers of events remaining after each selection step are summarized in 
Table [H We found 317 candidate events in the final data sample. Three types 



t3 



Additionally, energy-dependent cut criteria were also applied. 



of background events are considered. 

Table 1 

Summary of the number of remaining events and the cut signal efficiency for electron 

scattering and inverse beta decay after each selection. 



Cut 


Data 


Detection eff. 


Detection off. 






{u + e) 


(Pe + p) 


Energy range cut, 








Fiducial volume cut, 


2108 


0.99 


0.99 


Noise event rejection 








Spallation event cut 


1774 


0.91 


0.98 


Cherenkov opening angle cut 


759 


0.86 


0.93 


7-ray cut 


738 


0.83 


0.90 


Pion-like event cut 


705 


0.81 


0.88 


Multi-ring cut 


686 


0.81 


0.87 


Decay electron cut 


317 


0.81 


0.87 



The first dominant source of background events are decay electrons from invisi- 
ble muons induced by the atmospheric (anti-)muon neutrinos. The second such 
source are electrons generated by the atmospheric (anti-)electron neutrinos, 
while the third are the multiple nuclear de-excitation gamma-rays produced 
by NC interactions of the atmospheric neutrinos. 

The remaining spallation events are estimated to be 4.7±4.3 events in SK-I/III 
and 0.2 ± 1.7 events in SK-II using the likelihood distribution. 

The number of remaining solar neutrino events is estimated to be one. 

To extract signal events, an angular distribution with respect to the Sun di- 
rection is used. Figure [3] shows the cos6'sun distribution of the 317 candidates, 
where ^sun is the angle of the candidate event direction with respect to the 
expected neutrino direction calculated from the position of the Sun at the 
event time. To fit the distribution, we use a x^ defined in Eq.(4): 



20 



{A^obs.i - N^eAh) - N^pAh) - aA^bkg.i}^ 



X(Io,«)=2. ■ —2 —2 ' (4) 

j=l '^stat,* ~r '^sys,! 

where we use A'"obs,i as the observed number of events in the z-th angular bin. 
Nue,i(J-o) ^^^ Nup,i(Jo) are the expected signal events from electron scattering 
and inverse beta decay in the i-th bin, where Iq is the total flux of all flavors 
of neutrinos in units of cm~^s~^ (Fig- II])- The theoretical calculation of the 
angular distribution of the inverse beta decay is from Ref . [19] . The systematic 
errors used in the x^ function (cTsys,* = const.) are summarized in Table [2l 

-^bkg,i is the area-normalized background events for the i-th bin. We deter- 
mine the A''bkg,j values using the final data sample by applying an iterative 
procedure (described in Appendix) to remove a possible bias caused by the 
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Fig. 3. Angular distributions with respect to the expected neutrino direction from 
the Sun. The points with error bars show the data events. The soUd and dashed 
histograms represent Iq = Igo and Iq = 0, respectively, which are described in the 
text. 
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Fig. 4. Expected angular distributions of elastic scattering recoil electrons (left) and 
inverse beta decay positrons (right) assuming a neutrino flux of 3.0 x 10^ cm~^s~^ 
and the detection efficiency of the cuts. 

solar-correlated events. The factor a representing the background normaliza- 
tion is determined so that the x^ value is minimized for each Iq value. In the 
fit, we restrict a > 0. As a result of the fit, we obtain Iq = 102 ± 64 cm~^s~^, 
corresponding to a 1.6 a excess. 



We determine the 90% CL. upper limit flux Igo by using the following equa- 



Table 2 

Systematic errors (%) used in the ^ function. The systematic error on the cross 

section of inverse beta decay is from Ref. [20 j . 



tion: 



Error source 



SK-I SK-II SK-III 



Cross section 


0.4% 


0.4% 


0.4% 


Veutrino propagation 


5.0% 


5.0% 


5.0% 


Reduction efficiency 


2.3% 


3.8% 


2.1% 


Fiducial volume 


1.3% 


1.1% 


1.0% 


Livetime 


0.1% 


0.1% 


0.1% 



Total 



5.7% 6.4% 5.6% 



|^'%xp(-x'(Io)/2)dIo 
y^ exp(-x'(Io)/2)dIo 



(5) 



Figure [5] shows the exp (— x^(Io)/2) distribution as a function of Iq, where the 
normahzation is scaled so that the maximum value becomes one. We obtain 
Igo = 166.6 cm~^s~^. Figure |3] shows the cos^sun distribution superimposed 
with curves representing Iq = Igo and Iq = 0. We obtain a x^(Io = 0) value of 
20.02 compared to 19 degrees of freedom, which corresponds to a probability 
of 39.3% and is therefore consistent with a hypothesis of null signal in the 
final data sample. 
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Fig. 5. Likelihood function (oc exp (— x^(Io)/2)) as a function of monopole-induced 
neutrino flux (Iq). The dashed line represents Iq = Igo- The deviation from the Gaus- 
sian shape beyond lo = 183.4 cm~^s~^ is due to the restriction of a > 0. 



The monopole-catalyzed proton decay rate /„ in the Sun is given by, 



fv= \ nuVrcicrpp^Ad^x decays/s. 



(6) 



10 



where nu is the monopole number density, v^ei = /3reic the relative velocity 
between the monopole and the proton, a = cTo//3rei (^o ~ 1 iiib) the catalysis 
cross section [TUIIllf21] . pp the proton mass density, and Na Avogadro's num- 
ber (6.0 X 10^^). The space integral is taken over the interior of the Sun. Helium 
gives a negligible contribution [21j . If we assume the monopoles are accumu- 



lated in the center of the Sun where pp = 50 g/cm^ and (3re\ = 1.7 x 10~v^ 
we obtain fp = 1.7 x 10^(cro/l nib)NM, where Nm is the integrated number 
of monopoles in the region. The rate fp is given by fp = ATrdHo/S/fi,^, where 
(i = 1.5 X 10^^ cm is the distance between the Sun and the Earth. The factor 3 
takes into account the fact that three neutrinos are emitted per proton decay, 
/j,^ = /^+(1 — aT,+ ) is the fraction of z/g produced in a proton decay, where /^+ 
is the branching fraction of a proton decay into 7r+ -t- anything, and is about 
0.5 for some GUT models. The factor a.,^+ = 0.2 is the absorption probability 
of TT"*" at the center of the Sun [9J. Using the Igo value, we obtain an upper 
limit on Nm at the 90% C.L.: 



^" (i^) (fe) < '^' ' ''" '^'"'"^■■^■' 



(7) 



The monopole flux is calculated from the following equation by assuming 
monopole-antimonopole annihilation is negligible: 

Nm = 7ri?| [ 1 + (^ ] I e(MM, /3M)47rFMt0, 

where /3esc = 2 x 10~^ and /3m are the escape velocity and the monopole velocity 
in units of the speed of light, respectively, Rq is the solar radius (7.0 x 10^° cm) 
and ^0 the elapsed time after the birth of the Sun (4.6 x 10^ yr). e(MM, /3m) is 
the efficiency with which the Sun captures monopoles that strike its surface, 
which is dependent on the monopole mass Mm and velocity /3m. We estimate 
e(MM, (^m) based on the calculation of the energy loss of a monopole inside 
the Sun from Refs. [22123] ■ From ([7]) and ([SD, we obtain for /3m < 10"^ 




< 6.3 X 10-'^ -^ cm-^s-isr-^ (90% C.L.) (9) 



This limit is only valid for Mm < 10"*^^ GeV since e(MM,/3M) decreases as 
the monopole mass increases and becomes significantly less than 1 for higher 
monopole masses. 

Our result is presented in Fig. [6l Also shown are monopole fiux limits obtained 



* We take t;rei as the thermal velocity of proton in the center of the Sun, /3reic = 
{2kBTc/mp)^''^ ~ 1.7x 10~^ m/s, where /cb is the Boltzmann constant, Tc = 1.5 x 10^ 
K is the temperature at the center of the Sun, and nip = 0.938 GeV is the proton 
mass. 
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Fig. 6. 90% C.L. upper limits on the monopole flux as a function of monopole veloc- 
ity, /3m- Each flux limit is referenced in the text. Catalysis cross section ((To=lmb) 
and monopole mass (Mm = 10^^ GeV) are assumed for indirect searches. 



by various recent direct detection experiments: UCSD-II [21], Soudan-2 [25], 
Baksan [26], Kolar [27], Kamioka track etch [28], Mica [29] and the MACRO 
final result [30], as well as ones from indirect detection experiments: the Parker 
bound [H], limits from neutron star observations [21], and the Kamiokande 
experimental results [521153] . As seen in the figure, our new limit is many 



orders of magnitude more stringent than those of the direct experiments and 
is also two orders of magnitude lower than the result of the Kamiokande 
experiment [32][33] . which used a similar detection method. This is mainly 
because SK is a larger detector and the background events are subtracted 



12 



effectively using angular distributions. The present limit is lower than the one 
obtained from the X-ray excess of old neutron stars, 

Fm (^) < 3r X 10"^^ cm~2s"^sr"\ (10) 

Vlmb/ 

where Pu is 10"'^ and r is the ratio of the total luminosity to the photon 
luminosity of old neutron stars. The value of r varies from 1 to 10^, where 
r ~ 1 corresponds to the conventional neutron star equations of state while 



large values of r correspond to the possibility of a pion condensate ^^ I or quark 
matter core [21]. /Srei in old neutron stars is of order of 0.1 ~ 0.3. Therefore, 
our limit is better than the neutron star limit even in the most stringent 
case. Finally, we emphasize that this result also provides constraints on the 
parameter space of inflation models as suggested in [7,^8] . 



5 Conclusion 



A search for monopole-induced neutrinos from the Sun in the electron to- 
tal energy range from 19 to 55 MeV was carried out using a data sample 
of 2853 days' livetime and 22,500 metric tons fiducial volume. We found no 
signal events and set the 90% confidence level limit on the monopole flux: 
FM(o-o/lmb) < 6.3 x 10"^''(/3m/10~^)^ cm"^ s"^ sr"\ Our result provides the 
world's most stringent upper limit for /3m < 10^^. 
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^^ If there is a pion condensate in a neutron star, its total luminosity is dominated 
by neutrinos generated in the following reactions: tt~ + n -^ n+i~ + i>i and n+i~ — > 
n + 7r~ + f£, where £ = e or /i. 
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Appendix A Estimation of {A^bkg,j} 



Although the background shape in the angular distribution with respect to 
the direction of the Sun is expected to be almost isotropic, we tried rejecting 
directional biases including those caused by solar-correlated events remaining 
in the final data sample, which will distort the cos ^sun distribution. 

We define Nfin as the number of the final candidates in each SK phase and 
prepare event direction vectors, di, and solar direction vectors, Sj for the final 
sample {i,j = 1, ■ ■ ■ ,Nfin). Then, the angle with respect to the direction of 
the Sun of the i-th event is expressed as cos6'sun,i = di ■ Sj. 

First we fill the Nf}n(Nf}n — l)/2 possible combinations of d^ ■ Sj {i 7^ j) into a 
histogram hi, and fit hi with a polynomial fi{x),x G [—1, 1]. Also, we fit the 
actual true cos^sun distribution with another polynomial c{x). 

Then we fill the Nfin (Nfin — l)/2 combinations of di ■ Sj {i 7^ j) into another 
histogram /i2 with weights of Wi{di ■ Sj) = fi{di ■ Si)/c{di ■ Sj). 

Again we fit /i2 with another polynomial f2{x). 

Usually, fi{x),Wi and hi sufficiently converge after several iterations of the 
above procedure. In this way, we eliminate the bias from solar-correlated events 
in the limit that fi{x) describes the histogrammed distribution. 
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